Coherence scale of coupled Anderson impurities 



O 
(N 

> 

O 

:^ 

o 
en 






I 

o 
o 



en 
> 

in 

in 

o 
o 



% 



Lijun Zhu and Jian-Xin Zhu 

Theoretical Division and Center for Nonlinear Studies, 

Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 

(Dated: December 1, 2010) 

For two coupled Anderson impurities, two energy scales are present to characterize the evolution 
from local moment state of the impurities to either of the inter-impurity singlet or the Kondo singlet 
ground states. The high energy scale is found to deviate from the single-ion Kondo temperature and 
rather scales as Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction when it becomes dominant. 
We find that the scaling behavior and the associated physical properties of this scale are consistent 
with those of a coherence scale defined in heavy fermion systems. 
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I. INTRODUCTION 
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There has been much recent debate on the nature 
of a coherence scale for heavy fermion systems, which 
characterizes the evolution of /-electrons from the local- 
ized magnetic moments to itinerant quasiparticles [l|, [2] ■ 
In experiments this scale is manifested as the temper- 
ature T* above which the quasiparticle signatures van- 
ish, for example, the Drude peak in the optical conduc- 
tivity disappears, the magnetic entropy becomes i?ln2, 
and the resistivity (from magnetic ions) reaches a max- 
imum value. Traditionally, this is associated with the 
single-ion Kondo temperature Tk which characterizes the 
Kondo renormalization of /-electrons to quasiparticles by 
forming spin singlet resonance state with conduction elec- 
trons. It is not known a priori that such a single impurity 
picture is applicable to the lattice, where magnetic ions 
are coupled by intersite interactions. This can be exam- 
ined, for instance, by tuning the concentrations of mag- 
netic ions in a parent non-magnetic compound while the 
dilute and dense limits provide information on the single- 
ion Kondo behavior and the Kondo lattice behavior with 
coupled local moments, respectively. One of such studies 
is on Cea;Lai_j;CoIn5 and it is found that T* increases 
from the Ce-dilute limit with the increase of the concen- 
tration X and can become one order of magnitude larger 
in the dense limit [3]. An extensive analysis on a group 
of heavy fermion compounds displaying quantum critical 
properties [l| shows that T* for these materials are indeed 
bigger than Tr-. It is further discovered that T* scales 
as PqJx, where po is the density of states of conduction 
electrons at the Fermi energy and Jk is the onsite Kondo 
coupling between spins of magnetic ions and conduction 
electrons. This is the scaling form for the intersite spin 
exchange interaction, namely, Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interaction /, rather than the single-ion 
Kondo temperature Tk ~ (l/po)e"^/^'"'"'^'^. 

To resolve the discrepancy between T* and Tk there- 
fore relies on the understanding of the interplay between 
the intersite coupling and the onsite Kondo coupling. 
While the associated lattice models for heavy fermion 
systems are in general difficult to solve, the two-impurity 
Anderson model (or the equivalent two-impurity Kondo 



model) presents such a competition effect in an exactly 
solvable way. With the assistance of the cluster dynami- 
cal mean-field theory (DMFT), a self-consistently solved 
two-impurity model also provides a solution to the An- 
derson lattice model 4:]. In the two- impurity Anderson 
(or Kondo) model, an antiferromagnetic RKKY interac- 
tion between two local moments favors an inter-impurity 
spin singlet state which competes with the Kondo sin- 
glet state, which is similar to the same type of compe- 
tition in heavy fermion systems ^, 161] . It is known that 
for ferromagnetic RKKY interactions and a range of an- 
tiferromagnetic RKKY interaction, the evolution from 
the (high energy) local moment state to the (low energy) 
Kondo resonance state or the inter-impurity singlet state 
is through a two-stage process, characterized by two en- 
ergy scales. Previous studies on this model were focused 
on the low energy scale, which characterizes the quan- 
tum phase transition or crossover between the two types 
of ground states |6l-[l4|. The high energy scale, however, 
receives little attention hitherto, which is the focus of this 
study. While there are indications that it may be related 
to the single-ion Kondo temperature or RKKY interac- 
tion scale, there are no systematic studies on its scal- 
ing behavior and its associated physical properties, espe- 
cially, when RKKY interaction is much bigger than the 
single-ion Kondo temperature. The reason to carry out 
this study is that, we show, this high energy scale shares 
the same properties as the coherence scale defined for the 
heavy fermion systems. Following the experimental defi- 
nition, we denote this high energy scale as the coherence 
scale for coupled Anderson(Kondo) impurities 0). 

Besides the difference in the focus on the coherence 
scale, our study differs from previous ones in the follow- 
ing aspects. 1) We calculate the T = uniform and stag- 
gered spin susceptibilities in addition to the spectral func- 
tions (the imaginary parts of the single-particle Green's 
functions). We carry out an extended analysis on these 
dynamical quantities, for instance, a scaling analysis, and 
employ different quantities to identify and determine the 
characteristic energy scales. Therefore, we provide a de- 
tailed account for the physical properties associated with 
these energy scales. 2) We calculate these dynamical 
quantities with the recently-developed complete-Fock- 



space numerical renormalization group method [15| . This 
method was developed to overcome the spectral weight 
loss in the intermediate and high energy range from tra- 
ditional patching scheme in dynamical quantity calcu- 
lations [16i] . which becomes severe in the two- impurity 
problem with a larger eigen space. This approach is es- 
sential for determining the coherence scale which falls in 
this energy range. In particular, we show by a compari- 
son that the results from the patching scheme may lead 
to a misleading understanding of the coherence scale. To 
our best knowledge, this is also the first practice to gener- 
alize this method to multi-impurity models. 3) It is also 
important to choose a system for which the single-ion 
Kondo temperature Tk and the inter-impurity spin ex- 
change interaction / can be easily determined or directly 
given, which the identified characteristic energy scales 
can be compared to. We present results for two systems 
where two impurities are located far away from each other 
and at nearest neighbors on a three dimensional cubic 
lattice. For two impurities located far away, RKKY in- 
teraction, which is generated from the virtual exchange 
of conduction electrons, is found to vanish. Theoreti- 
cally, we can add a direct spin exchange term to sim- 
ulate RKKY interaction effect [y, 0). This system has 
been studied earlier, and we repeat the calculations on 
this system because it provides a reference system for the 
single-ion Kondo physics where Tk can be determined. 
We also provide new results for / 3> Tk- For two impu- 
rities sitting on nearest neighbors, a finite antiferromag- 
netic RKKY interaction between two impurities is found 
to be generated perturbatively in the order of poJfc, and 
is antiferromagnetic. As Tk and / are tuned in the same 
fashion as in realistic heavy fermion systems, this system 
provides a direct comparison to the experiments on the 
scaling behaviors of the characteristic energy scales. It 
also reveals an origin for a parity-splitting term which is 
inherent to the heavy fermion systems. 

From our numerical studies on the above two systems, 
we can in general identify two characteristic energy scales 
for two coupled Anderson impurities. While the ground 
state is found to be always a Fermi liquid fixed point, 
either a Kondo resonance state or the inter-impurity sin- 
glet state controlled by RKKY interaction, the low en- 
ergy scale Tl serves as the (local) Fermi liquid scale where 
Fermi liquid behaviors emerge. The high energy scale Th 
is found to be a spin fluctuation scale, where the imag- 
inary parts of the spin susceptibilities reach their maxi- 
mum values: the origin can be the Kondo spin-fiip scat- 
tering or inter-impurity singlet to triplet excitation. In 
either case, the single-particle (charge) excitation gains 
considerable weight. We further find that Th for both 
systems has the same behavior: it increases from Tk 
and then becomes / when I /Tk increases. This can 
be explained by whether the Kondo renormalization can 
reach Tk or is already cutoff by the inter-impurity sin- 
glet to triplet excitation gap. Our results indicate that 
the physical properties above Th are still determined by 
the single-ion Kondo physics, which allows us to make an 



argument that Th can be generalized to the lattice sys- 
tem. Th ~ I and its physical properties imply that it is 
consistent with the coherence scale T* defined for heavy 
fermion systems. For completeness, we also show the be- 
havior of Ti, which is the characteristic energy scale for 
the phase transition or crossover and is consistent with 
previous studies. Tj, has different behaviors in the two 
considered systems due to two different origins, the prox- 
imity to a degeneracy point of the two competing grounds 
states and a parity-splitted quasiparticle hopping term. 
The rest of the paper is arranged as follows. In SecHIl 
we introduce the two-impurity Anderson model and the 
numerical method we adopted. The results for the di- 
rect and generated RKKY interaction cases are shown 
in Sec. IIIIl where two energy scales are identified and 
their properties are characterized. We focus our discus- 
sions on the high energy scale, in Sec. IIVI where we show 
that the scaling behavior of this scale and its associated 
physical properties are consistent with those of the co- 
herence scale identified in the heavy fermion systems. In 
Appendix, we provide some technical details on how to 
evaluate Tk and RKKY interaction. 



II. MODEL AND METHOD 

The Hamiltonian for the two- impurity Anderson model 
can be written as 

+ E ^//-/- + E Unf.nnf,i + /Si • S2 , (1) 

ia i 

where i sums over two impurity sites. This describes two 
interacting local orbitals fia (Anderson impurities) em- 
bedded in a non-interacting conduction electron medium 
Ckrr with the system size Nc and in hybridization with 
Ckff with the strength Vk at each impurity site r^. e/ is 
the local orbital energy level and U is the onsite Coulomb 
interaction for the impurities. / is a direct spin exchange 
interaction between two impurity spins. In reality, such 
term, known as RKKY interaction, is perturbatively gen- 
erated through a virtual process exchanging conduction 
electrons. In the single occupancy limit of each orbital, 
this model can be mapped into a two-impurity Kondo 
model (some details are presented in Appendix). 

By taking the even or odd parity combinations of the 
local orbitals fp=(e,o),cr = (/ict±/2o-)/'\/2, the fluctuations 
due to conduction electrons can be represented by two 
separate baths with the hybridization functions 
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In our calculation we assume Vif^ — V and a three dimen- 
sional tight-binding dispersion et = ~{D /i) X]i=i *^os fc^a 
for conduction electrons with D the half-bandwidth and 



a the lattice constant. We consider two cases with impu- 
rities located far away [case (i)] and at nearest neighbor 
lattice sites [case (ii)], where jri — r2| = cx) and a, re- 
spectively. Te.oi^) for these two cases can therefore be 
determined with Eq. ([2]). Subsequently we can determine 
the Kondo coupling Jk and the generated RKKY inter- 
action / while the calculation details are presented in Ap- 
pendix [^ We find that the hybridization functions for 
case (i) and case (ii) can be written as Te^oi^) = To and 
Tf.^oi'-^) = TqCITw), respectively, where Fq — npoV^. We 
take U = 2D and e/ — —U/2 in both cases. The Kondo 
coupling constant is found to be pqJk — STq/^ttU). The 
generated RKKY interaction vanishes in case (i), indi- 
cating that it is effectively a single-ion Kondo problem, 
from which we can determine Tk. Theoretically, we add 
a direct spin exchange term /Si -82 to simulate RKKY in- 
teraction. In case (ii), an antiferromagnetic RKKY inter- 
action is perturbatively generated, which is determined 
to be / « 0.20pqJ]^ (see Appendix [X| . Therefore we 
no longer need to add a direct spin exchange term. The 
ratio I/Tk can be tuned by varying the Kondo coupling 
since / and Tr- have different dependencies on Jk- In our 
calculation, this is achieved by varying V as Jk ~ V'^/U. 

This impurity problem is amenable to the numerical 
renormalization group (NRG) method 16] , which pre- 
scribes a non-perturbative procedure to capture the low 
energy properties. It discretizes the hybridization func- 
tion into a chain of electron orbitals with the energy 
decreasing in logarithmic scale. In the two-impurity 
problem, as protected by the parity symmetry, Fg and 
To can be discretized separately. While only the head 
site of the chain couples to the impurities, we can solve 
the original model by an iterative diagonalization to in- 
corporate gradually more low energy sites. In prac- 
tice, as the eigen space increases with more sites in- 
cluded in each iteration, only certain number of low en- 
ergy states are kept for the next iteration. From these 
NRG eigenstates, we can calculate dynamical quantities, 
UAsiuj) = -iJ^dte"^*{[A{t),B{0)]±), in particular, the 
Green's function Gfpa with A = B'^ = /po-, the uniform 
(x«) and staggered {xa) impurity spin susceptibilities, 
with A = B = {Si, + S2,)/V2 and {Si, - S2,)/V2, re- 
spectively. Traditionally, these quantities are calculated 
from the kept states from all iterations. However, kept 
states from different iterations are not necessarily orthog- 
onal to each other. A patching scheme is adopted to avoid 
the overlap contributions, for instance, in Ref. [7|. This 
is to choose a weight function for the overlap spectrum 
between different iterations. The arbitrariness of the 
weight function leads to the inaccuracy of these dynam- 
ical quantities. The recently developed complete-Fock- 
space (CFS) method 15| overcomes this problem, which 
we adopt in our calculation. It instead calculates dynam- 
ical quantities from all the discarded states (including all 
states from the final iteration), which is found to form 
an orthogonal Fock space conserving the total spectral 
density. Although the calculation time almost doubles 
with an additional backward iteration to determine the 



reduced density matrix, we find that the improvement 
in dynamical quantities is significant, especially in the 
energy range near and above Tk, which is of interest to 
our study. We will show a comparison for these meth- 
ods in Fig. [6] While we typically choose a discretization 
parameter A = 2 and keep 4000 lowest energy states in 
each iteration, we find that by choosing a smaller A or 
keeping more states, the dynamical quantities have little 
change in the high energy range. 



III. EMERGENCE OF TWO ENERGY SCALES 

A. Case (i): two impurities located far aw^ay from 
each other with a direct spin exchange interaction 

We first consider case (i), two impurities located far 
away from each other which can be represented by 
Te{oj) = Toioj) = Fq. Although this case has been stud- 
ied earlier [7[, we re-examine it here for two considera- 
tions. Firstly, as shown below, it provides a reference 
system to the single-ion Kondo physics when the direct 
spin exchange term is not added. Tk can be determined. 
Secondly, it is a rare system exhibiting a quantum criti- 
cal point with the added direct spin exchange term, with 
or without particle-hole symmetry (this is related to the 
vanishing parity-splitting potential scattering term). 

The results of the spectral function Af{uj) = 
— ImG'/(a;)/7r, imaginary parts of the uniform (x") and 
staggered (x") spin susceptibilities are shown in Fig. [T] 
We choose a particle-hole symmetric case e/ — — C//2 ~ 
-D, and Fe,o = 0.0457rD for jwj < D. Here D = 1 
serves as the energy unit. The unit for spin susceptibili- 
ties {gpsY is also set to 1, where g and p,B are the Lande 
factor and the Bohr magneton, respectively. 

In this case, without the direct spin exchange inter- 
action /, we find that Xu{<-^) — Xa('^), i-e., there is no 
correlation between two impurities {SizS2z) = 0. There- 
fore, the spin dynamics, which controls the low energy 
properties, is the same as in the single- impurity Ander- 
son model. In finite U case, the Anderson model dif- 
fers from its Kondo counterpart in additional features 
in high energies, i.e., the incoherent peaks of A/(w) at 
a; « £/ and e/ + U corresponding to the free-orbital 
fixed point. In general, there exists an additional energy 
scale characterizing the evolution from the free-orbital 
fixed point to (or close to) the local moment fixed-point, 
where charge fluctuations freeze and Af{uj) is small. In 
the Kondo limit, [/—>■€», these two fixed points can be 
thought of being pushed to the infinite energy. The low 
energy physics is controlled by the Kondo renormaliza- 
tion from the local moment fixed-point to the strongly- 
coupling Kondo fixed point, and can be characterized by 
a single energy scale, the single-ion Kondo temperature 
Tk- Below Tk, Af has a resonance peak centered at 
the Fermi energy. [v4/(0) = l/(7rFo) by Friedel's sum 
rule. Our result has a few percent deviation which is 
due to the truncation of eigen space in NRG. In practice. 
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FIG. 1: (color online) Spectral functions (a), imaginary parts 
of the spin susceptibilities (b) as functions of energy for vari- 
ous values of / in case (i) , two impurities located far away but 
coupled with a direct spin exchange interaction /. x" ^nd Xa 
in (b) are represented by dotted and solid lines, respectively. 
Here Te = To = 0.0457rD and e/ = -(7/2 = -D. While only 
the a; > range is shown, it is found that Af{—Lj) = Af{ij) 
and x"{~'^) ~ ~x"{^)- ^f is also symmetric for different 
parities and spins in this case. Two energy scales Th and Tl 
(see definitions in the text) are illustrated for I/D = 0.002 
and 0.01. 



we find that the improvement for this low energy sum 
rule can be achieved by keeping more states in NRG it- 
erations and/or adopting a two-particle Green's function 
method 17]]. Tk is also a spin-fluctuation scale, indi- 
cating the spin fluctuations, originated from the Kondo 
spin-flip scatterings, reach the maximal strength. As the 
a; = Kondo fixed point is a Fermi Iqiuid fixed point 
with scattering phase shift de,o = 7''/2, (local) Fermi 
liquid behaviors, such as Af{uj) — Af{0) ^ —uj^ and 
x"a('^) ~ ^^ emerge when |w| < Tp (Tp is smaller 
than Tk but is not an independent scale as it is pro- 
portional to Tk)- We determine Tk from the real part of 
the spin susceptibilities, x'u ai'^ = 0) — l/(4?i<:), which 
can be determined from the imaginary parts with the 
Kramers- Kronig relation or the Korringa relation [18| : 
lim^^o[x::,a('^)H = Cu,a = 27T[x'u,aiO)?- We find that 
;^;^_^(0) = 2.5 X 10^, and 2.7 x 10^ by these two relations, 
and we determine pqTk ~ 1.0 x lO^'^. 

Upon adding an antiferromagnetic spin exchange in- 
teraction /, the low energy properties have changed: 
the spectral function is reduced from an energy scale 
around Tk, then then either increases again to the full 
weight l/(7rro) or decreases to [a pseudogap state with 
Af{uj) ^ uj"^, verified by a log- log plot] at a lower en- 
ergy scale. The abrupt spectral weight change at the 



Fermi energy is a signature of the two-impurity quantum 
phase transition [a,13|: the ground state changes from a 
Kondo resonance states with two impurity spins forming 
Kondo singlets with the conduction electrons, to an inter- 
impurity spin singlet states with two impurity spins form- 
ing singlets by themselves and being decoupled from the 
conduction electrons. The critical values is found to be 
Ic ~ 2.ZTk which is also consistent with previous stud- 
ies [g, 0| . Apparently, the single- ion Kondo picture with 
a single scale, Tr-, no longer applies and two energy scales 
are needed to explain the dynamical quantities. On both 
sides of the critical point, Fermi liquid behaviors persist 
at low energies, such as 'x'u,a{'^) = Cu,a^- We therefore 
associate the low energy scale T^ as the Fermion liquid 
temperature. But we need to pay attention that Tj, alone 
cannot differentiate the Kondo resonance state and the 
inter-impurity singlet state, both of which are Fermi liq- 
uid fixed points but have different spectral weights or the 
scattering phase shift at the Fermi energy. The high en- 
ergy scale Th is associated with a spin-fluctation scale, 
where x'u ai'^) reach the maximal strength (for I close to 
Ic, where Xal*^) becomes relatively a constant, Th is the 
high energy edge). 

According to this observation, we define Th as the 
peak position of x"(w) and T^ as the onset energy when 
Xu ai'^) becomes linear in w. The former can be easier 
determined from the numerical data. The errors come 
from systematic ones such as numerical discretization and 
the broadening of the delta function to a log-Gaussian 
form. The latter, however, has arbitrariness because 
the crossover to the Fermi liquid behaviors spans a wide 
energy range. Therefore, we also provide characteriza- 
tions from l/[4Xu.a(0)] ^''^^ Cu,a- In the single-ion case, 
l/[4x^,a(0)] = Tk while C„,a = Cq - l/T^. When / 
is close to Ic, we find that the constant piece in x"(w) 
remains relatively unchanged, implying CaT^ — CqTk- 
In this region, TKCo/Ca provides a more rehable method 
to determine Tl- We plot the energy scales determined 
from these methods in Fig. [21 

From Fig. [21 we learn that Tl and TKCo/Ca have 
similar behaviors (the magnitude difference is similar to 
the difference in Tp and Tk in the single- ion case). As 
TkCq/Cu is determined directly from numerical data, it 
shows better asymptotic behaviors. Near Ic, we find that 
it can be fitted as {I — Ic)^- This indicates a uniformly 
vanishing energy scale, which characterizes the continu- 
ous quantum phase transition. The exponent 2 is also 
consistent with previous studies [3, 0] ■ But away from 
Ic, the exponent is smaller or bigger than 2 for / — >■ 
and I :^ Ic respectively. In the latter case, it is rather 
fitted as {I-Icf/^ within a range oi I - Ic > 0.2/^. This 
fractional exponent may not be universal, as we can also 
fit 1/Ca -- l/Cn ~ /* for / » Ic- Another check for 
the low energy scale is from a scaling analysis, by rescal- 
ing the energy with respect to 1/Ca- We carry on such 
an analysis for Xa(w), which is shown in Fig. [31 Once 
the energy is rescaled with TKCo/Ca, the low energy 
part of Xai^) f^-Us in the same universal curve within 
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FIG. 2: (color online) (a) Characteristic energy scales as func- 
tions of I/Tk- The energy scales determined from 1/Xu,a(0), 
and Cu,a ~ liin;j_»o x",a('^)/'^ a-re also shown, (b) shows 
the low energy scale as a function of (7 — Ic)/Ic- Two 
dotted lines are fitting lines T/Tk = 2.0(///c — 1)" and 
T/Tk = 4.0(///c - 1)'*''^. (c) shows amplified region of the 
high energy scale. Two dotted lines are the guide line T = I 
(black) and a fitting line T/Tk = 0.8{I/Tk)^ for 1/C„. 



\6I/Ic\ < 20%, which is associated with the universal- 
ity of the two-impurity quantum critical point. Indeed, 
TjiCo/C'a provides a faithful representation of the low en- 
ergy scale Ti . On the other hand, the Korringa relation 
is violated for finite Is. While 1/Xa(0) has the same trend 
as Tl, it rather vanishes logarithmically ^ 1/ In |/ — Id- 
Th and l/[4x^(0)] have similar behaviors: they in- 
crease uniformly when / increases. TkCq/Cu follows the 
same trend, but increases more rapidly when I > Ic- 
Clearly, Th increases from Tk and becomes the RKKY 
scale / when / ^ Tk- 



B. Case (ii): two impurities sitting on nearest 
neighboring sites with generated RKKY interaction 

We proceed to study case (ii), two impurities sit- 
ting on nearest neighbors which can be represented by 
Te^oi^) = ro(lT'^)- The results are shown in Fig.|4l We 
observe similar behaviors, where small (large) V cases 
can be compared with large (small) I /Tk cases in case 
(i). The differences are, Af for even and odd parties 
split due to different Fe^ol'^), and A/(0) smoothly varies 
rather than a jump. We can also identify two similar 
energy scales in Aj and x" a and plot them as functions 
of V in Fig. Sfc). As learned from case (i), both Th 




FIG. 3: (color online) Scaling analysis on the staggered spin 
susceptibility. x"('^) as functions of a; are shown for different 
Is, as in Fig. [TJb). Here ui is scaled with TKCo/Ca to the 
scaling behavior. 



and Tl deviate from Tk when RKKY interaction is fi- 
nite, it is not possible to determine Tk from the dy- 
namical quantities in this case. Instead, we calculate 
Tk from the corresponding single impurity case with 
re.o(w) = Trpo^^ as in case (i) with the same V . We find 
that it can be fitted by the standard expression pqTk = 
0.5exp[-(l/poJK) + (l/2)ln(poJK) + 1.58(poJif)'] where 
PoJk = 8poV'^/U. RKKY interaction is determined to 
be antiferromagnetic / « 0.20poJj^ (see Appendix [X)) . 
As V increases, both / and Tk increase monotonically 
but their ratio decreases. This provides a direct simula- 
tion to realistic heavy fermion materials: both follow the 
Doniach's phase diagram [cf. Fig. Hl^c)] [19| . 

Since Tk and Fq are varying as functions of V in this 
case, it is helpful to replot our results in scaled forms. 
For instance, we can rescale energies with Tk and Af 
with l/(7rFo) for each value of V. This allows a direct 
comparison to case (i) where Tk is fixed. The rescaled 
spectral functions and energy scales are plotted in Fig. \5\ 
For the spectral functions, it is interesting to observe that 
the high energy parts follow the universal Kondo curve, 
i.e., neglecting the incoherent peaks, this curve describes 
the Kondo renormalization from a local moment state 
at infinite energy to the low energy Kondo fixed point. 
However, this curve is rather associated with a particle- 
hole asymmetric case, or with a finite potential scattering 
term. The potential scattering terms for the even and 
odd parities have the same magnitude and are opposite 
in sign. The potential scattering term in the single-ion 
model is known not to change the universality of the 
Kondo fixed point: it simply shifts the resonance peak 
position and contributes an additional phase shift tt/2 + 
Sp. Therefore, Af{0) is reduced from l/(7rFo). The V = 
0.4D case with / ^ Tk presents such an example. When 
V decreases, I/Tk increases, and the low energy physics 
begins to be affected by RKKY interaction. Similar to 
the particle-hole asymmetric two-impurity model studied 
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FIG. 4; (color online) Spectral functions (a), imaginary 
parts of the spin susceptibilities (b) as functions of energy 
for various values of hybridization constant V in case (ii), 
two impurities sitting on nearest neighbors with generated 
RKKY interaction. The solid (dashed) lines represent the 
the even (odd) channels for Af, and Xa iXu) for spin sus- 
ceptibilities. With Te{—ijj) = Foiijj), Ae,o{i^) = Ao,e{—Uj) 

is satisfied, (c) shows different energy scales as functions 
of V (in unit of D). Tk is determined from calculations 
by assuming Te,o = To, and is fitted (blue dashed line) by 
PoTk = 0.5exp['--(l/poJK) + (1/2) ln(/9o Jx) + l.58{poJKf]. 
The black dashed line is the RKKY scale I = 0.20po Jx. 



before: the system shows a crossover from the Kondo 
resonance state to the inter-impurity spin singlet state, 
rather than a continuous transition. Tj, remains finite 
near / w 2.2Tk (at V/D w 0.35). The divergence in 
Xa(0) is also absent. Nevertheless, we find that Th has 
the same behavior as in case (i): it uniformly increases 
from Tk to / when I /Tk increases. 



C. Origin of two energy scales 

The emergence of two energy scales indicates the ef- 
fect of a finite RKKY interaction on the Kondo renor- 
malization. While the triplet configuration of the two 
impurity spins, as an 5 = 1 impurity, can undergo the 
spin-fiip scattering to itself in favor of Kondo resonance, 
the singlet configuration only couples to the conduction 
electron by being firstly excited to the triplet configura- 
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FIG. 5: (color online) The rescaled spectral functions (a) and 
energy scales (b) in case (ii), impurities sitting on nearest 
neighbors with generated RKKY interaction. These are the 
same plots as in Fig. [3] (a) and (c). After a rescaling, they 
can be compared with Fig. [TJa) and Fig. (2] 



figurations of the two impurity spins, causing these two 
configurations to renormalize differently. This can be ev- 
idenced from the difference in the uniform and staggered 
spin susceptibilities. When RKKY interaction is antifer- 
romagnetic [20] but small, though the singlet configura- 
tion is favored, the Kondo scattering between the singlet 
and triplet configurations, on the energy scale of Tk-, is 
able to overcome the singlet-to-triplet gap and the Kondo 
resonance state persists at low energies. When RKKY in- 
teraction is much larger than T/f, the Kondo resonance is 
suppressed at low energies, leading to a state with gapped 
quasiparticle excitations, denoted as the inter-impurity 
singlet state. 

While both the Kondo resonance state and the inter- 
impurity singlet state belong to the Fermi liquid fixed 
points, Ti characterizes the energy difference between 
these two states. Above T^, the system evolves into 
an intermediate state where the inter-impurity singlet to 
triplet excitation and the Kondo scattering are on equal 
footing and the physical properties deviate from FL be- 
haviors. This intermediate state is also the quantum crit- 
ical state in case (i) when T^ vanishes a,t I = Ic- From our 
results, we learn that the spectral weight for the quantum 
critical state is exactly the half of that for the Kondo reso- 
nance state and the staggered spin susceptibility diverges 
logarithmically. We find that these properties are con- 
sistent with the Majorana fermion picture proposed by 
the Bosonization approach (lOl . Il3 | . This picture dictates 
that the quantum critical state is a partially screened 
Kondo state with only half of the impurity degrees of free- 
dom as a Majorana fermion forming a resonance model 
with the extended electron degrees of freedom. This can 
also be evidenced from the thermodynamics calculated 
from NRG ^1] that the entropy (per impurity) for this 
state is (ln2V2. While Tl ~ {I-Ic^/Tk in case (i) indi- 
cates the proximity to the degeneracy point between the 
Kondo resonance state and inter-impurity singlet state, 
the finite Tl in case (ii) is due to another effect. In case 
(ii), the particle-hole symmetry is actually broken in each 
channel due to the presence of a parity-splitting potential 



scattering term V-{c\^Cea — c\„Cocr)- It is shown [3, [i3| 
that such term will generate a coupling between the other 
half of the impurity degrees of freedom (the other flavor of 
Majorana fermion) and the extended degrees of freedom. 
Therefore, there is always a finite quasiparticle spectral 
weight at the Fermi energy. The FL temperature is also 
finite Tl - TK{poV-f [1,111, [12. We notice that such 
a parity-splitting term can also be generated by other 
forms of particle-hole asymmetry, such as a hybridiza- 
tion (hopping) term tf{fl^f2a + h.c.) 7], the asymmetry 
in energy dependence of Kondo coupling jl3| or a regu- 
lar potential scattering term [iJl- Compared with these 
studies, case (ii) we studied here points out an origin for 
the parity-splitting term which is inherent to the heavy 
fermion systems, and cannot be tuned away. 



IV. THE RELATION BETWEEN Th AND THE 
COHERENCE SCALE 



While the properties of T^ are the focus of previous 
studies on this model [6l-[l4| and our results are in agree- 
ment with these studies, we focus our discussions on scal- 
ing behavior and the properties of Th which has received 
little attention. 

As shown in Figs.[ljc) and[SJb), Th has the same scal- 
ing behavior in the two cases studied above: Th ~ Tk 
(or /) for / ^ Tk (or ^ Tk) and a value in-between 
when / and Tk are comparable. When / ^ Tk, the 
Kondo effect dominates. A small RKKY interaction in 
this limit only enhances slightly the staggered spin fiuc- 
tuations and affects the physical properties in the en- 
ergy range (Tk — I,Tk) while the Kondo renormaliza- 
tion above Tk remains unaffected. Therefore, Th ~ Tk 
in this limit. Tk lies in the crossover region between the 
local moment and the strong coupling fixed points in the 
single-ion Kondo renormalization. Correspondingly, the 
imaginary parts of the spin susceptibilities change behav- 
iors at Tk from ^ l/w to ^ w, i.e, the spin fiuctuations 
reach the maximum strength associated with the spin-flip 
Kondo scattering. The quasiparticle spectral weight also 
increases rapidly at Tk- Indeed, Th (being Tk) serves 
as the onset scale for the formation of coherent quasipar- 
ticles. We notice that quasiparticles at this energy scale 
do not necessarily exhibit Fermi liquid behaviors, which 
rather develop at a lower energy scale. When / 3> Tk, 
RKKY interaction suppresses the Kondo scattering from 
the Fermi energy up to the scale /. The Kondo renormal- 
ization is cutoff by I before reaching Tk- The imaginary 
parts of the uniform and staggered spin susceptibilities 
reach the maximum values at /, but have a hump struc- 
ture instead of the smooth crossover between the local 
moment to FL behaviors as in the Kondo resonance state. 
While this hump is due to the singlet to triplet spin fluc- 
tuations between two impurity spins, it is broadened by 
the Kondo scattering. In other words, impurities are still 
weakly coupled to the conduction electrons to form quasi- 
particles, i.e., the coupling is Jk rather than JK.eff -^ oo 



at the Kondo flxed point. The spectral weight is still fi- 
nite and reaches the maximum value at / though it is 
small. Th is determined by / in this limit, which can be 
understood as a cutoff effect on the Kondo renormaliza- 
tion. When / is comparable with Tk, the Kondo scatter- 
ing and inter-impurity singlet to triplet excitation have 
comparable strengths. Th as the location for the max- 
imum spin fluctuations lies in between Tk and / and 
increases when I/Tk increases. In all these cases, we 
flnd that Th serves as the scale where the physical prop- 
erties begin to be strongly influenced by the formation 
of quasiparticles. Therefore, we can associate Th to the 
coherence scale for two coupled impurities. 

We notice that the properties of Th are in agreement 
with those of the coherence T* in heavy fermion systems. 
For instance, the changes in optical conductivities, the 
tunneling spectroscopy, and the Hall effect are related to 
the change of the quasiparticle spectral weight while the 
magnetic responses and the scattering rate are affected 
by the strong spin fluctuations at this scale. As observed 
from the two-impurity model, the single-ion Kondo tem- 
perature Tk loses its manifestation in physical properties 
when / > Tk- This is also true for heavy fermion sys- 
tems. In experiments, Tk is determined by the dilute 
magnetic ions limit, as in Ref. [3|. Indeed, it is in gen- 
eral one order of magnitude smaller than T* [l|. This 
implies that most heavy fermion compounds fall in the 
/ > Tk regime with magnetic orders developing at low 
temperatures. In this case, we flnd that Th is indeed de- 
termined by RKKY interaction, which is consistent with 
the experimental analysis [l|. In this case, we notice 
that the quasiparticles at the scale / are different from 
Landau quasiparticles displaying Fermi liquid behaviors: 
they rather are underscreened Kondo excitations with 
non-Fermi liquid behaviors. Nevertheless, the dynam- 
ical measurements determining the coherence scale in- 
deed are not necessary in the Fermi liquid regime. Be- 
sides, Fermi liquid behaviors develop at a lower energy 
scale. We flnd that the Fermi hquid parameters, such 
as Cu.a — linioj^oix" a('^)/'^]; are also determined by /. 
In this sense, the coherence scale indeed appears in the 
Landau parameters, which are manifested also in static 
property measurements. This provides an explicit expla- 
nation to the scaling behavior of T* as RKKY interaction 
rather than the single-ion Kondo temperature, and it can 
be interpreted as a coherence scale. 

Then the question is whether Th, the coherence scale 
for two Anderson impurities, can be generalized to the 
heavy fermion lattice systems. This relies on an ex- 
plicit calculation on the lattice models, for instance, solv- 
ing the periodic Anderson lattice model through a self- 
consistent two-impurity Anderson model with the cluster 
DMFT 4]. In this approach, one chooses the neighboring 
two sites of the lattice as a cluster and treat the lattice 
as the repetitions of this cluster. The even and odd par- 
ity states therefore correspond to the degrees of freedom 
at momentum (0,0,0) and (7r,7r, tt), respectively, for a 
three dimensional lattice. Once self-consistently solved. 
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FIG. 6: (color online) Comparisons between the spectral func- 
tions obtained by CFS-NRG method (solid lines) and the 
patching scheme (dotted lines) (please refer to Sec. [IT] for 
details of these two methods). The parameters taken are 
the same as in case (i), or Fig. [T] For both methods, the 
NRG iteration part is the same. A weight function w{x) — x 
is commonly chosen for the patching scheme. The spectral 
functions obtained by CFS-NRG method have high weights 
in the whole energy range. We explicitly verify the sum rule 
J dLjAf{iu) = 1, and find the deviation in practice is less than 
0.1%. 



the two impurities have the same properties as the two 
sites in the lattice incorporating two-site correlations. Al- 
though the self-consistency procedure hasn't been done 
explicitly, we argue that the generalization of the coher- 
ence scale can be valid. The most important result from 
our calculations is that the high energy properties above 
Th are governed by the single impurity physics, i.e., the 
intersite coupling only modifies the low energy proper- 
ties. This can be observed from the spectral functions 
[Cf. Fig. (l)a and Fig. 3(a)], which follow the universal 
single-ion Kondo curve at high energies but their low en- 
ergy weights are suppressed systematically from the scale 
Th, i.e., Th being RKKY interaction, is a cutoff scale on 
the single- ion Kondo renormalization. This statement, 
however, cannot be made from the patching scheme cal- 
culations. In Fig. [SI we show the spectral functions ob- 
tained from both CFS-NRG method and the patching 
scheme. While these methods may give the same posi- 
tions for the characteristic energy scales, as the NRG iter- 
ation part is the same, the spectral function obtained by 
the patching scheme may lead to a different understand- 
ing on Th for / 3> Tk- We observe enhanced quasipar- 
ticle excitations at w sa / compared with the single-ion 
case. This is due to the choice of the weight function be- 
tween the overlapping energy range, which, commonly, 
favors the states in later iterations (lower energies) and 
exaggerates its weight. However, this may indicate that, 
due to RKKY interaction J, considerable spectral weight 
transfers from the incoherent part to the scale / to form 
coherent excitations, rather the cutoff effect indicated by 



the CFS-NRG result. It is from the spectral functions 
that we determine the new Weiss fields to recalculate the 
two- impurity problem. As the high energy part above Th 
is determined by the single-ion Kondo physics, it remains 
unchanged in the self-consistency procedure. The low 
energy part, however, is expected to be modified in the 
self-consistency procedure. But as / is the cutoff scale, 
it does not involve any spectral weight transfer from the 
incoherent part, or more strictly, from part with energy 
higher than /. Therefore, Th remains as the high energy 
scale for the lattice. A caveat for this argument is not 
to violate the spectral sum rule. As seen from our cal- 
culation, the spectral weight of the low energy part only 
consists a negligible portion of the total spectral weight, 
for instance, valid for Ce compounds. We notice that 
this argument may not hold for cases where the weight 
transfer between the incoherent part and coherent part 
is significant. 



V. CONCLUSION 

In summary, we have studied the two-impurity An- 
derson model with the complete-Fock space numerical 
renormalization group method. We find that due to the 
competition between the Kondo effect and RKKY inter- 
action which is the spin exchange interaction between two 
impurities, there are in general two energy scales present. 
The high energy scale characterizes the quasiparticle for- 
mations and the low energy scale specifies when these 
quasiparticles develop Fermi liquid behaviors. In be- 
tween, there exists an intermediate state with non-Fermi 
liquid behaviors where both the Kondo scattering and 
the excitations between the singlet and triplet configu- 
rations of the two impurity spins play equally important 
roles. 

The focus of our study is on the high energy scale which 
we associate with the coherence scale in heavy fermion 
systems. In the two-impurity model, it marks the scale 
with maximum spin fiuctuations. While the spin fiuctua- 
tions can be due to the Kondo scattering with conduction 
electrons or the singlet to triplet excitations between the 
impurity spins, their relative strength is tuned by the 
ratio between RKKY interaction and the Kondo tem- 
perature. When the Kondo effect wins, the high energy 
scale is determined by the Kondo temperature, where 
quasiparticles begin to gain considerable weight. When 
RKKY interaction dominates, the high energy scale is de- 
termined by RKKY interaction, which can be understood 
as a cutoff on the Kondo renormalization by a spin ex- 
citation gap. The quasiparticle spectral weight therefore 
reaches the maximum value at this scale. In both cases, 
the high energy scale characterizes when the quasipar- 
ticles gain considerable weight to influence the physical 
properties. This is consistent with the definition of the 
coherence scale in heavy fermion systems. The most im- 
portant result is that the physical properties above the 
high energy scales is determined by the single-ion Kondo 



physics and is expected to remain unchanged in a clus- 
ter DMFT procedure. This enables us to associate the 
high energy scale directly to the coherence scale of heavy 
fermion systems. Besides sharing the same physical prop- 
erties, they also have a same scaling behavior. When 
RKKY interaction is much bigger than the single-ion 
Kondo temperature, the coherence scale is determined 
by RKKY interaction which is consistent with the exper- 
imental analysis |l| . This provides an explicit theoretical 
description for such an energy scale. 



The low energy scale, as the Fermi liquid tempera- 
ture for quasiparticle, is expected to change in the clus- 
ter DMFT self-consistent procedure. This prevents a di- 
rect connection between the two- impurity quantum phase 
transition to the magnetic quantum phase transition in 
heavy fermion systems. However, we expect both ori- 
gins for the low energy scale in the two-impurity model 
will play important roles in determining the low energy 
scale (s) in a lattice model. The degeneracy point between 
magnetic excitations and Kondo scattering also exists in 
a lattice, for which Doniach's criterion remains the same, 
i.e., Ic/Tk is of order unity. If the two- impurity critical 
dynamics drives a magnetic transition for the lattice, it 
will be interesting to check the existence of Majorana 
fermion: many of its properties are in accordance with 
the local quantum critical behavior [23|. On the other 
hand, the two-impurity quantum critical point is unsta- 
ble against certain forms of particle-hole asymmetry. For 
instance, the parity-splitting potential scattering term is 
always finite in case (ii) of our study, which is inherent 
to the heavy fermion lattices. In a lattice, this can be 
related to asymmetry between degrees of freedom at mo- 
mentum (0,0,0) and (tt, tt, tt). We have shown that it not 
only generates a finite Fermi liquid temperature but also 
generates a finite spectral weight at the chemical poten- 
tial [cf. FiglUa)]. The parity-splitting term also relates 
to a finite Rel]i2(w — 0), which corresponds to a kinetic 
energy term for /-electrons [Ij]. This might help to sta- 
bilize a finite effective Kondo scale. Incidentally, as the 
potential scattering term scales as the Kondo coupling, 
the associated (low) energy scale then also has the same 
scaling relation as RKKY interaction. 



We expect that the two-impurity Kondo physics will 
also help to understand the Mott transition and uncon- 
ventional superconductivity in other strongly correlated 
systems. The Kondo physics has provided an elegant ex- 
planation to the Mott transition in the single-site DMFT 
approach to the Hubbard model [24]. However, the inter- 
site spin exchange interaction, similar to RKKY interac- 
tion, is absent in the single-site approach. We learn from 
the two-impurity model that such an intersite interaction 
competes with the local Kondo dynamics to modify the 
low energy properties. Indeed a multi-site cluster DMFT 
calculation shows different properties the metallic and 
insulating phases near the transition point [25|. 
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Appendix A: Determination of the hybridization 
functions and the interaction parameters 

We explain some of the calculation details for the hy- 
bridization functions and the interaction parameters such 
as the Kondo coupling Jk and RKKY interaction I . 
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FIG. 7: (color online) The hybridization functions for the 
even and odd parity channels for two Anderson impurities 
in a three dimensional cubic lattice, (a) is for two impurities 
located far away from each other (here ri — r2 = AQax) and (b) 
is for two impurities sitting on nearest neighbor sites (ri — r2 = 
ax) . Here V = 0.3D. 

The hybridizations functions can be calculated from 
Eq. ^ by the given hybridization constant Vk, the con- 
duction electron dispersion et, and the impurity loca- 
tions. We choose Vk = V , e\s_ — ~{D/3) J2i=i cos kia. 
For two impurities located far away from each other and 
sitting on nearest neighbors, we determine Te^oij^) nu- 
merically and the results are shown in Fig. [71 In Fig.[71[a), 
we find that Te^ui) = Toiuj) = irV'^pioj), where p{uj) is 
the density of state function for the dispersion Ck. Near 
the Fermi energy, re_o(w) are constants. In Fig.[71^b), we 
find that Te.oi'-^) ~ 1 T (^ near the Fermi energy. For 
simplicity, we take their low energy forms to the whole 
band. The difference is above the Van-Hove singularity 
at |a;| « 0.31? which doesn't affect the low energy physics. 

The Kondo coupling constant Jk can be obtained from 
a canonical transformation from the Anderson model to 
the Kondo model. In the two-impurity Kondo model, as 
the Kondo scattering can involve electrons from the same 
or different parity channels, there are different types of 
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Kondo exchange interaction 

V(w,w')/^ = - [rpiLo)rp,iuj')f^cioj,Lu') , (Al) 

TT 

where {p,p') = e, o and 

C(a;,a;') = — ^—— 

uj-ef uj-Ef- Uf 

1 1 

(A2) 



f'/ 



^' - e/ - Uf 



In comparison to the notations in Ref. 14], Je = Jee, 
Jo = Joo and J„i = (Jeo + Joe)/^. With e = -U/2 
and |a;|,|a;'| < C/, C{u},u}') « 8/C/. The Kondo couphng 



is taken as the value at the Fermi energy Jk = (J^ + 
J^ + 27^)1/2/^2 = 8ro/(7rpot^)- However, the energy 
dependence of the Kondo couphng is crucial to determine 
the generated RKKY interaction, which is 



2 



dw 



dbj' 



Ijj' — UJ 

(A3) 
If we take the hybridization functions to be constants, we 
find that / = — 21n2po(«/e — Jo)^- It is then always fer- 
romagnetic except when J^ = Jo, i-e., the two impurities 
located far away. When re,o(w) = ro(l=Fw), we evaluate 
/ « 0.20po'^I-, i-e., antiferromagnetic as expected. 
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